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We observe two-body loss of 3 Po 87 Sr atoms trapped in a one-dimensional optical lattice. We 
measure loss rate coefficients for atomic samples between 1 and 6 that are prepared either in a 
single nuclear-spin-sublevel or with equal populations in two sublevels. The measured temperature 
and nuclear spin preparation dependence of rate coefficients agree well with calculations and reveal 
that rate coefficients for distinguishable atoms are only slightly enhanced over those of indistinguish- 
able atoms. We further observe a suppression of excitation and losses during interrogation of the 
1 So- 3 -Po transition as density increases and Rabi frequency decreases, which suggests the presence 
of strong interactions in our dynamically driven many-body system. 



Ultra-narrow optical transitions in alkaline-earth 
atoms are the foundation for state-of-the-art optical lat- 
tice clocks [THS] and recent proposals for quantum infor- 
mation science (6j[7]- All of these applications will benefit 
from achieving long atom-light coherence with large num- 
bers of atoms, although at large densities inelastic and 
elastic collision processes become increasingly important. 
For ensembles of indistinguishable fermions, s-wave col- 
lisions are forbidden by the Pauli exclusion principle and 
higher order partial wave collisions are suppressed when 
the thermal energy of the atoms is much less than the 
energetic barrier for interactions. However, inhomoge- 
neous excitation can populate multi-particle electronic 
states that are not symmetric with respect to exchange, 
thereby allowing s-wave collisions (Sj [9] . Furthermore, p- 
wave collisions, though suppressed at low temperatures, 
are not forbidden and evidence of p-wave collisions was 
recently observed in an optical lattice of fermionic Yb 
atoms [TO] , 

In this Letter, we demonstrate that nuclear-spin- 
polarized, fermionic 87 Sr atoms in the 3 Po state experi- 
ence inelastic p-wave collisions that cause two-body pop- 
ulation decay in a one-dimensional (ID) optical lattice. 
Atoms prepared in an incoherent mixture of two nuclear 
spin sublevels experience slightly greater loss than atoms 
polarized to a single nuclear spin sublevel. We also ob- 
serve evidence of strong interactions during spectroscopy 
of the 1 S'o- 3 Po transition as we reduce the Rabi frequency, 
fl. Typically, a system is termed strongly interacting 
if the thermally averaged mean interaction energy per 
particle, U, dominates over all other energy scales. In 
our dynamically driven system, the finite temperature 
only affects the dynamics through a slight perturbation 
to based on the vibrational mode of an atom [5J 
Since the only other relevant energy scale is Ml, our sys- 
tem becomes strongly interacting when U/h ^> fi. Pre- 
viously, this was accomplished by confining atoms in a 
two-dimensional (2D) optical lattice, which increases U 
but decreases lattice site occupancy to mainly one or two 
atoms [TTJ [T2] . In this work, we demonstrate the pres- 



ence of strong interactions by decreasing fl in a ID lat- 
tice, where we can achieve an average of over 20 atoms 
per site. 

At a fixed density, strong interactions are signaled by 
an inhibition of losses and excitation with decreasing Q. 
As f2 decreases, the number of inelastic collisions dur- 
ing excitation increases for constant pulse area and an 
increase in losses is naively expected. Similarly, for suf- 
ficiently small we measure that both losses and exci- 
tation fraction are suppressed by increasing density be- 
yond a critical value, contrary to the expectation that 
loss will increase with density. Similar inhibition mech- 
anisms have been observed in other strongly interacting 
systems [TBHTo] . 

Both strong elastic and inelastic interactions could be 
responsible for the observed suppression in our experi- 
ment. Although losses are significant, a mean field den- 
sity matrix formalism that only includes inelastic pro- 
cesses underestimates the observed suppression. This to- 
gether with the observation of a density dependent asym- 
metric broadening of the lineshape suggest that elastic 
interactions are playing an important role in the dynam- 
ics. 

In our experiment, 87 Sr atoms are cooled to about 2 /iK 
in a magneto-optical trap (MOT) based on the 1 So- 3 -Pi 
transition [17] . The atoms are loaded into a vertically 
oriented, ID optical lattice that is overlapped with the 
MOT. The distribution of atoms across lattice sites is de- 
termined by the vertical extent of the MOT which is ap- 
proximately Gaussian with a standard deviation, a = 30 
/im. At the largest achieved density, we estimate that 
over 8000 atoms are distributed over about 400 sites with 
an average occupation number of 23. The size of the atom 
cloud in a lattice site is calculated from the temperature 
and lattice trap frequencies, measured using Doppler and 
sideband spectroscopy [9]. Since the extent of the 1 Sq- 
3 Pi MOT is small compared to the 3 mm Rayleigh range 
of our lattice beam, we assume identical trapping poten- 
tials for all sites. 

Atoms can be sideband and Doppler cooled (or heated) 
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FIG. 1. (Color Online) a) Atom number as a function of 
lattice hold time for atoms prepared in the g state (black 
triangles) and in the e state (green squares). Fits to the decay 
curves are calculated for the case of exponential decay (dashed 
lines). The solid line is a fit to the e atom decay using a sum 
over decays in single lattice sites [Eqn. (2)]. b) Comparison 
of e atom decay in the presence of an equal population of g 
atoms (purple circles) to that of a pure population of e atoms 
(green squares). Here, temperature is 3.5 fiK. 



on the 1 5'o- 3 Pi, F = 11/2 transition. Simultaneously, 
atoms are either pumped to the nip = +9/2 state (po- 
larized) using circularly polarized light resonant with the 
1 Sq- 3 Pi, F = 9/2 transition, or pumped to an incoher- 
ent mixture of spin states with equal populations in the 
mp = ±9/2 states (dual spin state) using linearly po- 
larized light on the 1 iSb- 3 Pi, F = 7/2 transition. The 
quantization axis is defined by a bias magnetic field par- 
allel to lattice polarization, which is strong enough to 
prevent nuclear spin depolarization. We quantify atom 
number in both 1 5'o (g) and 3 Pq (e) states by detecting 
fluorescence on the strong 1 So- 1 Pi dipole allowed transi- 
tion transition at 461 nm both before and after e atoms 
are re-pumped to g. To measure atom loss from the e 
state, atoms are excited on the g-e transition prior to 
lattice hold time with a resonant 7r-pulse from an ultra- 
stable laser co-propagating with the lattice beam and 
7r-polarized. The remaining g atoms are removed with a 
5 ms pulse of 461 nm light. 

Measured loss from the g state is well represented by 
exponential decay with a lifetime of 7 - 8 s, consistent 
with loss due to collisions with background gas. In con- 
trast, we measure a rapid, density-dependent loss from e 
that is inconsistent with a simple exponential decay law. 
Fig. [I] a) shows measured atom number as a function of 
lattice hold time for polarized samples of g and e atoms 
under similar temperature and trapping conditions. The 
additional loss from e results from inelastic e-e collisions. 
Fig. [I] b) compares e atom decay with and without an 
equal number of g atoms present. Agreement between 
the two curves in Fig. [I] b) limits inelastic g-e collisions 
to below the sensitivity of our experiment and thus, we 



neglect them. The decay of g atoms in the 50/50 mixture 
is likewise unperturbed by the presence of e atoms. 

To quantify loss from e, we adopt a model that includes 
both one and two-body losses. The atomic density in e, 
n e , is described by 



K, 



(1) 



Here, T is the one-body loss rate due to collisions with 
background gas, and K ee is the two-body loss rate co- 
efficient. As in Refs. [HI [19], spatial integration of the 
solution to Eqn. ([I]) yields an expression for the atom 
number in a single lattice site as a function of time: 



N(t) 



N exp(-Tt) 



1 + [N Q K ee /(TT^Tw^w z )}[l - exp(-rt)] ' 



(2) 



Here, No is the initial atom number in site and w z (w r ) is 
the 1/e 2 radius of the atom cloud in the strongly (weakly) 
confined direction(s). 

For polarized atoms, loss occurs dominantly from in- 
elastic p-wave collisions at microkelvin temperatures. 
This gives that K ee = K p nd {T), where K™ d (T) is the 
loss rate coefficient due to inelastic p-wave collisions be- 
tween indistinguishable e atoms, which depends on tem- 
perature, T . For dual spin state atoms, loss can occur 
from intra-spin-state odd partial wave collisions and from 
inter-spin-state collisions which can be any partial wave. 
Keeping only s and p-wave contributions, the decay of a 
single spin state, a, in the presence of another spin state, 
P, can be written as, 



n a = -Tn a - K;r(T)ni - (K? st + K dist {T))n a n ti , (3) 

where n a (rip) is the density of spin state a (j3), and _ft. dlst 
is the loss rate coefficient due to ^-wave inelastic collisions 
between distinguishable e atoms. In the case where n a = 
np = 1/2 n e the differential equation for the total e state 



1/2 K 



dist 



3/4K p nd (T) 



density gives that K ee 
Kl nd = 2Kf ist [201 mi- 

To compare with experiment, we extract K ee from the 
decay in total atom number using a fit to a sum of single 
site decays [Eqn. ([2])] based on the estimated distribu- 
tion of atoms. The one-body decay rate is set to 1/7.6 
s _1 , the measured value from g atom decay. Rate coef- 
ficients for dual spin state and polarized atoms at differ- 
ent temperatures are shown in Fig. [2] One might expect 
to measure much greater rate coefficients for dual spin 
state atoms compared to polarized atoms since p-wave 
collisions should be suppressed at microkelvin tempera- 
tures. Yet, measured loss coefficients for dual spin state 
atoms are only slightly larger than for polarized atoms 
at equivalent temperatures. 

To understand this result, we perform a time- 
independent quantum calculation, similar to that in 
Refs. [221 [23], using a single scattering channel, and a 
short-range boundary condition at an interatomic sep- 
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FIG. 2. (Color Online) Measured two-body loss rate coeffi- 
cients for polarized and dual spin state atoms as a function of 
temperature. The vertical error bars are calculated from an 
uncertainty in atom cloud size associated with a temperature 
uncertainty of 0.5 fjK, a conservative estimation to account 
for measurement uncertainty and experimental drifts. Solid 
lines are calculated values for loss rate coefficients (see text). 



aration, R — R Q , described by two parameters. A first 
parameter, 8, represents an accumulated phase-shift from 
R = to R = Rq due to an unknown atom-atom short- 
range potential of Sr2. A second parameter, p\ s , rep- 
resents the probability of two atoms to be lost when 
they encounter at R — R . These atomic losses are due 
to couplings that can take place between different elec- 
tronic potential energy curves of the Sr2 complex at small 
R [24H27] . The release of kinetic energy associated with 
changes in electronic configurations results in trap loss. 
The long-range interaction potential between two 87 Sr 
atoms is given by an attractive van der Waals electronic 
potential. We choose an isotropic C% van der Waals co- 
efficient of 5260 a.u. (1 a.u. = 1 Ehag, Eh is the Hartree 
energy, arj is the Bohr radius) for the e-e interaction |28j . 

The logarithmic-derivative of the scattering wavefunc- 
tion is computed for each R after giving an initial value 
at R = Rq which is defined in terms of the two pa- 
rameters 8 and pi s . Using asymptotic boundary con- 
ditions at large interatomic distances, we obtain cross 
sections for a wide range of collision energies. Thermal- 
ized loss rate coefficients are calculated by averaging the 
cross section over a Maxwell-Boltzmann distribution of 
the relative velocities in three dimensions. Using values 
of 8 = 0.9 ir, pis = 0.4 and Rq = 30 ao as initial boundary 
conditions, we were able to simultaneously determine: (i) 
K™ d = Tx(4± 2) 10~ 6 cm 3 s" 1 K" 1 (ii) and Kf ist = 



for 87 Sr atoms, as well as re- 



(ra 2 x 10~ n cm^s" 1 ) and (iv) the 
10" 



(1.4 ± 0.8) 10- 11 c 
produce (iii) -?Q nd 

elastic cross section (« 7 x 10~ 12 cm 2 ) of indistinguish- 
able bosonic 88 Sr e atoms, from previous experimental 
studies [19J . From a collisional point of view, the loss 
rate probability at short range pi s = 0.4 indicates that 
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FIG. 3. (Color Online) Measured Rabi lineshapes of g to 
e interrogation for different experimental conditions. Black 
circles show the fraction of population in the e state and green 
triangles show the total number of both g and e atoms. In 
panels a) and b) f2 = 2n x 2.5 Hz, corresponding to a 200 
ms 7r-pulse. All experimental conditions in these two panels 
are equivalent except that in panel a) the atom number is 
much greater than in panel b). Panels c) and d) are similarly 
related but for Q — 2n x 1 Hz, corresponding to a 500 ms 
7r-pulse. 



the Sr-Sr system deviates significantly from a high-lossy 
universal system (pi s = 1) [221 El] where s-wave collisions 
are generally two orders of magnitude higher than p-wave 
collisions [H] at these typical temperatures. 

To investigate the effect of inelastic collisions on Rabi 
spectroscopy of the g to e clock transition, we interro- 
gated the transition with an ultra-stable laser capable of 
resolving sub-Hz spectral features [30] with varying atom 
numbers and values of tt. In these measurements, lattice 
trapped atoms are cooled to 2 /iK and polarized. For each 
fi, the probe time is held constant and the pulse area is 
controlled using laser power to our best approximation of 
a 7r-pulse by optimizing excitation on resonance in a low 
density condition. We vary the number of atoms loaded 
into the lattice by adjusting the loading rate into a first 
stage MOT based on the 1 S'o- 1 i : i transition. The line- 
shape for each experimental condition is a superposition 
of ten individual scans. To eliminate the effect of laser 
drift on the measured lineshapes, we cancel the residual 
drift of our ultra-stable laser (about 80 mHz/s) to less 
than 5 mHz/s. Furthermore, we alternate the direction 
in which we scan across resonance. 

Fig. [3] demonstrates the dramatic effect that atomic 
interactions have on Rabi spectroscopy of 87 Sr. One dis- 
tinct feature in the measured lineshape for Q = 2?t x 2.5 
Hz at the largest achieved density [see panel a)] is a side- 
band at negative detuning. Similar features have been 
observed in a tightly confined 2D lattice [12]. In that 
work, the observed sidebands were attributed to inho- 
mogeneous excitation in doubly occupied sites which al- 
lows transfer to electronic states that are antisymmetric 
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with respect to exchange and therefore separated from 
nonintcracting electronic states by the energy of s-wave 
interactions. The main spectral features in Ref. [12] were 
always dominated by lattice sites with only one atom and 
the possible role of p-wave interactions in doubly occu- 
pied sites was not considered. The key difference in this 
work is that most atoms are in multiply occupied lattice 
sites. Therefore, a systematic evaluation of the role s 
and p-wave elastic collisions play in spectroscopy of the 
g-e transition will be necessary before the origin of these 
spectral features is conclusively determined. 

The lineshapes in Fig. [3] clearly demonstrate the sup- 
pression mechanisms present in our system. The differ- 
ence in excitation fraction between vertically adjacent 
panels demonstrates the suppression of excitation frac- 
tion with increasing density. Furthermore, the difference 
between a) and c) demonstrates the suppression of both 
excitation fraction and loss for similar densities as 11 de- 
creases and confirms that loss is also suppressed as den- 
sity increases for sufficiently low 11 since loss does not 
increase in c) compared to in d). 

Although a quantitative analysis of these lineshapes 
has the potential to illuminate the collisional processes 
at work in our system, their asymmetric features indi- 
cate the presence of strong interactions and will require 
a true many-body approach. Nonetheless, we have mod- 
eled the peak excitation fraction and atom loss measured 
at different densities and 11 using a mean field density 
matrix formalism, similar to that in Ref. |18j . In this 
model, the density matrix, p, within a lattice site evolves 
as 

p = ~[h, p ]+h(p), (4) 

where H is the atom-light Hamiltonian and 7Z(p) is a 
relaxation matrix which describes population and coher- 
ence decay. In the rotating wave approximation these 
can be written as 

f = (n/2 T)' np)n = -r Pll 
n(p) 12 = K( P y 21 

= - {K ee p 22 /2 + T/2 + L + K dep (p n + pxd) Pia 

Tt{p) 2 2 = - (r + K ee p 22 ) P22, 

(5) 

where A is the detuning from resonance and K& cv is a 
phenomenological rate coefficient that accounts for elas- 
tic dephasing collisions. L is set to 0.05 x 17 to account for 
excitation inhomogeneity arising from the dependance of 
11 on trap oscillator level [S]. 

We determine the peak e population and peak atom 
loss from each binned lineshape and express both as a 
fraction of the measured atom number. The excitation 
fraction at a particular detuning is calculated using the 
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FIG. 4. (Color Online) a) Peak excitation (black circles) and 
peak atom loss (green triangles) vs. Rabi frequency. Atom 
loss is expressed as a fraction of the total atom number in 
the absence of loss (off resonance), b) Peak excitation and 
peak atom loss vs. atom number for a fixed Rabi frequency 
of 2n x 2.5 Hz. Solid lines show calculations using the den- 
sity matrix model described in the text with Ka ep inversely 
proportional to the Rabi frequency. Dashed lines indicate the 
same calculation with the only difference being Kdep is set to 
zero. 



total atom number at that detuning, while the atom loss 
compares the atom number at the detuning with great- 
est loss to the mean atom number away from resonance. 
In Fig. [3J these quantities are plotted versus 11 keeping 
atom number constant, a), as well as versus atom num- 
ber for a fixed 17 of 2tt x 2.5 Hz, b). The dashed lines are 
calculations of the plotted quantities using the density 
matrix model above with Kd cp set to zero and K ee set 
to the maximum value allowed by our measurement at 
these experimental conditions. Solid lines correspond to 
another calculation with iQep set to C x (2n x 5 Hz /SI) 
with C = 3.5 x 10 -12 cm 3 /s. Both K ee and Kd cp are 
necessary for experiment-theory agreement. 

The variation of -Kdcp with SI is necessary to achieve 
agreement with the experimental data and is motivated 
by the analytic solution to Eqn. Q with K ee and L set 
to zero. In the limit where 0/(n e Xdep) R the solu- 
tion predicts a suppression of excitation fraction propor- 
tional to Sl/i^dcp = Q 2 /C, a form similar to the predicted 
behavior for suppression mechanisms in other strongly 
interacting systems [TTJ [T3]. This suggests that elastic 
interactions are playing an important role in the dynam- 
ics. Nevertheless the addition of this phenomenologi- 
cal parameter produces lineshapes that are broader than 
those experimentally measured and thus a full many- 
body treatment is needed for understanding the role 
played by inelastic and elastic processes. 
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Supplemental material 



j P Loss Data 



We quantify atom number by detecting fluorescence 
on the strong 1 S a — 1 P\ dipole allowed transition at 461 
nm. For each experimental condition, we make a series 
of 30-50 repeated measurements of the atom number at 
various lattice hold times and calculate the mean and 
standard error for each series. Preceding each series of 
measurements at a particular hold time we perform a se- 
ries of measurements with zero lattice hold time. In order 
to reduce the effect of long term drifts in atom number, 
the mean value for a particular hold time is expressed as 
a fraction of the mean value with no additional lattice 
holding. These fractional values are scaled by the mean 
value of all zero hold time measurements and the uncer- 
tainty of this value is calculated from the uncertainty of 
the fractional value and the standard deviation of the 
zero hold time measurements. 

In order to extract the loss rate coefficient from our 
measured loss curves we first estimate the distribution of 
atoms in the lattice. This distribution is determined by 
the vertical extent of our second stage MOT based on the 
So — 3 P\ transition. The MOT cloud is approximately 
gaussian with a standard deviation of approximately 30 
/an. We assume that the number of atoms loaded into 
a lattice site obeys a Poissonian probability distribution 
with a mean value determined by the assumed distribu- 
tion of atoms, 



X(x) 



V2 



Ntot 

^ ex P 



7T<7 Z 



V25- 



(A.6) 



Here, N to t is the measured number of atoms, a is the 
standard deviation of the gaussian distribution in units of 
lattice sites [a = a/(X/2) for A the lattice wavelength and 
a the standard deviation in units of length] , and x is the 
number of lattice sites from MOT center. By summing 
the Poissonian distribution in atom number and lattice 
sites over all sites 5a from center we obtain a value for the 
number of lattice sites containing N atoms for each value 
of N from 1 to 100. Even for the largest achieved samples 
of atoms the average number of lattice sites with more 
than 60 atoms is much less than 1. The experimental 
data is fit to a sum of decay equations of the form of Eqn. 
(2) from the main text, ranging in initial atom number 
from 1 to 100. Each term in this sum is weighted by the 
number of atoms and the average number of lattice sites 
with that number of atoms. 

The spatial extent of the atom clouds are assumed to 
be gaussian and identical across lattice sites since all oc- 
cupied sites are well within the 3 mm Raleigh range of our 



lattice beam and we assume a uniform temperature. Us- 
ing the measured temperature and trap frequencies, we 
estimate the 1 /e 2 radius of the cloud in the i direction to 
be 



Uh 



X y/2fa) + 1, 



(A.7) 



where m is the mass of the atom, Vi is the trap frequency 
in the i direction, and (n,) is the average vibrational 
quantum number in the i direction, given by 



(e hv >' kBT -l) \ 



(A.8) 



where ks is the Boltzmann constant. 

The extracted value of K ee is determined by a fit to 
the data with each data point weighted by the inverse 
square of its uncertainty. We estimate that our tempera- 
ture measurements are accurate to ±0.5 /iK and the final 
uncertainty in K ee is dominated by the variation of the 
fitted value when the temperature is changed by ±0.5 
/iK, leading to variations in the Wi's. 



Lineshape Data 

To acquire a single measurement of the Rabi lineshape, 
we step the frequency of an ultra-stable laser across the 
resonance frequency of the 1 S'o - 3 Po transition. At each 
frequency point, we measure the number of atoms re- 
maining in the Sq (g) state by detecting fluorescence 
on the 1 Sq — 1 P\ transition while the atoms are being 
excited on resonance by a probe laser. The heating as- 
sociated with this measurement removes all (g) atoms 
from the lattice. Atoms in the 3 Po (e) state are then re- 
pumped into g using 679 nm light resonant on the 3 P 
— Si transition and 707 nm light resonant on the 3 P2 - 
3 Si transition. These re-pumping lasers are broadened 
to about 1 GHz to account for hyperfine splitting of the 
3 P~2 and 3 Si states. We then measure the number of 
the re-pumped atoms. An additional probe laser pulse is 
applied to measure the background light level. 

From these three quantities, we calculate the number 
of atoms as N tot = N g +N e — 2*Nb, where Nb is the back- 
ground level and N e (N g ) is the number of e (g) atoms. 
The excitation fraction is then N e /N to t- To reduce sta- 
tistical fluctuations in the measured lineshapes, we align 
multiple scans by a Lorentzian fit to the lineshapes and 
superpose them. The resulting lineshape is binned into 1 
Hz or 0.5 Hz bins and the excitation fraction of each bin 
is determined from the mean of the points within that 
bin. The uncertainty is determined from the standard 
error of the mean. 



